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Motivated by recent experiments on rotating Bose-Einstein condensates, we investigate a rotating, 
polarized Fermi gas trapped in an anharmonic potential. We apply a semiclassical expansion of the 
density of states in order to determine how the thermodynamic properties depend on the rotation 
frequency. The accuracy of the semiclassical approximation is tested and shown to be sufficient 
for describing typical experiments. At zero temperature, rotating the gas above a given frequency 
Qdo leads to a 'donut '-shaped cloud which is analogous to the hole found in two-dimensional Bose- 
Einstein condensates. The free expansion of the gas after suddenly turning off the trap is considered 
and characterized by the time and rotation frequency dependence of the aspect ratio. Temperature 
effects are also taken into account and both low- and high-temperature expansions are presented for 
the relevant thermodynamical quantities. In the high-temperature regime a virial theorem approach 
is used to study the delicate interplay between rotation and anharmonicity. 

PACS numbers: 51.30.-hi, 31.15.xg, 67.85.Lm 



I. INTRODUCTION 



In a remarkable series of experiments in the nineties, sophisticated cooling techniques were used to expose the 
quantum nature of ultracold dilute gases. For bosons, the realization of Bose-Einstein condensation (BEC) in 1995 
[1-3] revealed the importance of the underlying aggregating statistics for the existence of a broken symmetry phase 
in weakly interacting systems. For fermions, the statistics has just the opposite effect. Due to the Pauli exclusion 
principle, low energy collisions of fermionic particles are strongly suppressed, worsening the thermalization between 
the atoms and, therefore, making the cooling techniques more involved. Nevertheless, a few years later another 
milestone experiment was realized in which a quantum degenerate Fermi gas was observed [4]. Since then, research 
with both types of degenerate quantum gases has experienced a strong growth, and different trends have been pursued, 
each showing interesting and relevant connections to other fields of physics such as atomic, molecular and condensed 
matter physics. 

Rotating quantum gases constitute one of these fast growing areas. In the case of bosons in a harmonic trap, 
intensive investigations have been carried out, both theoretically [5-11] and experimentally [12-14]. If the rotation 
frequency Q is of the order of (but smaller than) the radial trapping frequency cc;_l, i.e. Q < cc;_l, the rotation is 
called slow. In this regime, intrinsic properties of rotating systems like the temperature dependence of the moment 
of inertia [15], rotation induced excitations [16], and the formation of vortices [17, 18] have been considered. If one 
rotates faster, vortex arrays are created [19] where small- amplitude oscillations, so called Tkachenko modes, have been 
observed [20] and theoretically investigated [21-24]. In addition, fast rotating {Vt ^ uj^) BECs have been linked to 
fascinating many-body phenomena such as the integer quantum Hall effect [25, 26] or its fractional version [27], and 
are suspected to display quantum phase transitions to non superfluid strongly correlated states [28, 29]. 

In the fast rotation regime, experiments in harmonic traps are limited by the loss of confinement due to centrifugal 
forces when the rotation frequency approaches the axial harmonic trapping frequency. To overcome this problem, A. 
L. Fetter suggested adding a quartic term to the trap potential [30]. The successful experimental implementation 
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of such a quartic term was carried out in the group of J. Dahbard in Paris [12]. In their set-up a laser is used 
to stir the quarticahy trapped condensates, making possible the exploration of many properties of the system such 
as, for instance, the vortex distribution, the optical density [12], and the shape oscillations [13]. The experimental 
realization triggered further theoretical research focusing on Bose-Einstein condensates (BEC). For instance, their 
collective oscillations [31], t her mo dynamical [32] and dynamical [33] properties have been considered in the presence 
of the quartic term. 

Thus, the question arises, what new features come about through the interplay between rotation and fermionic 
statistics. In a harmonic trap, some properties have already been investigated. For instance, the moment of inertia 
and its relation to quadrupole oscillations, in a framework valid for bosons and fermions [34], and the collective 
excitations [35] have been studied. In the superfluid phase, formation of vortices [36] has been observed and a 
rot at ion- induced phase separation has been proposed [37]. In addition, the influence of the Landau levels on the 
density profiles of two-dimensional polarized fermions in the context of the integer quantum Hall effect [38] and the 
experimental feasibility of fractional quantum Hall states [39] have also been theoretically considered. 

In the case of an anharmonic trap, apart from studying the breathing mode in the Bardeen- Cooper- Schrieffer 
(BCS)-BEC crossover regime [40], fermions have not been considered and most of the interesting aspects of the phy- 
sics of rotating systems remain restricted to pure harmonic confinement. At low temperatures, polarized atomic Fermi 
gases are essentially non-interacting due to Pauli blocking and the fact that p-wave interactions are negligible in such 
systems. However, due to the Fermi pressure, particles tend to avoid the center of the trap, so that, effectively, a 
statistical interaction arises even in the case of an ideal gas, which bears some analogies with repulsively interacting 
bosons [41]. 

In this work we are concerned with rotating ideal Fermi gases of spin polarized particles rotating with frequency 
Q. We consider the same trap configuration as in the Paris experiment, which will be described in detail in Section 
II. In Section HI we outline a general approach to formally justify the semiclassical approximation for calculating the 
density of states in anharmonic traps. Furthermore, we apply this technique to the Paris trap and see that a treatment 
within the semiclassical approximation is justified. In Section IV we consider ground-state properties of such a system 
in a similar semiclassical approach. We start by evaluating the Fermi energy for arbitrary rotation, then we obtain 
the particle density, the momentum distribution, and free expansion time-of-flight absorption pictures. In Section V 
temperature effects are considered and we derive expressions for the grand-canonical thermodynamic properties in 
terms of incomplete Fermi functions. In addition we analyze analytic low- and high-temperature expansions. For the 
latter, a virial theorem approach is also presented. Section VI discusses the main results of this work emphasizing the 
effects of rotation and anharmonicity in comparison to a non-rotating harmonic trap and presents an outlook for the 



In the following we consider a trap potential which is anharmonic and consists of harmonic axially symmetric 
magnetic trapping supplemented by an additional Gaussian trapping term in the xy-plane. The Gaussian laser beam 
is approximated by its fourth order Taylor series 



where Uq and w are the laser intensity and width, respectively. The radial distance to the origin is denoted by 
r± = The first term in (1) is just a shift in the energy scale and will be neglected in the following. The 

second order term will be absorbed into the harmonic trapping term. The last term is the one responsible for the 
quartic confinement and is characterized by the anharmonicity constant k = 2Uq/w^. It is important to notice that 
the harmonic trapping in the 2;-direction is not affected by the Gaussian laser potential and is characterized by its 
harmonic oscillator frequency cj^. 

An issue in the field of rotating quantum gases is that the rotation frequency, with which the condensate effec- 
tively rotates, might be different from the one of the stirring laser, due to changes in the shape dependent rotation 
properties of the gas such as the moment of inertia [14]. We neglect these effects and assume that a stirring laser 
can adiabatically add angular momentum to the particles, so that the gas rotates with frequency yielding an 
additional centrifugal harmonic term to the trap potential in the co-rotating frame [7, 8]. Under these considerations, 
the rotating anharmonic potential takes the form [32] 



field. 



II. ANHARMONIC TRAP POTENTIAL 
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FIG. 1: (Color online) The anharmonic potential (2) as a function of radial position in the plane ^ = 0, for the parameters 
given at Table I at both sub- and supercritical rotations. For the latter, a minimum at / appears due to the centrifugal 
nature of the harmonic term. 

TABLE I: Experimental values of the parameters of the anharmonic trap considered in this work. 
^ (Hz) (Hz) A k (Jm-^) n N 
64.8(3) 11.0 6 2.6(3) x 10"^^ 1.9 3 x 10^ 



where = and = ^JhjMuj^ are the axial harmonic energy and length scales, respectively, A = uo^/uOz denotes 
the anisotropy of the trap, n = ka^/cz represents the anharmonicity of the trap, and M is the mass of the atomic 
species under consideration. The rotation is taken into account with the parameter r] = 1 — Q'^ /uj'j_. Adopting 
the values given for the Paris experiment [12], the resulting rotating trap is then described by the radial harmonic 
trapping frequency u;± = 27r x 64.8(3) Hz, the axial harmonic trapping frequency cOz = 27r x 11.0 Hz, the force constant 
k = 2.6(3) X 10~^^ Jm~^ of the quartic trapping term, the anisotropy A ^ 6, and a tunable rotation frequency Q 
induced by the stirring laser. 

In the Paris experiment, a BEC was obtained with the bosonic isotope ^^Rb, for which the dimensionless anhar- 
monicity takes the value hi ~ 0.4. Since fermionic species like ^^K and ^Li have been sympathetic cooled by using 
the bosonic atoms ^^Rb [42, 43] and ^^Na [44], respectively, we consider them likely candidates for experiments with 
ultracold degenerate Fermi gases in the Paris trap. For definiteness, we evaluate all results of this paper for the 
parameters of the trap of the Paris experiment and the ^^K species, for which we have k ~ 1.9 (see Table I). 

Let us discuss the trap potential (2) and some of its peculiarities concerning rotation. For the moment, we restrict 
ourselves to the x?/-plane. As depicted in Fig. 1, for slow rotations, i.e. Q < u;±{r] > 0), the minimum of the potential 
is located at the origin. Rotating faster one eventually reaches the point where the harmonic confinement vanishes 
{Q = co±^r] = 0) and the trap is purely quartic. Hereafter this will be called critical rotation. If the rotation is made 
even faster, one achieves the s upercritica l regime, where the harmonic confinement becomes centrifugal and the trap 
minimum is moved to r± = ^J —r]y?a\l k and, with respect to the origin, has a depth given by 

l\E = \^rfez/4.K., (3) 

which naturally characterizes the anharmonic trap. The lowest available potential energy is, thus, —AE^ in the 
supercritical case. 

III. SEMICLASSICAL APPROXIMATION FOR DENSITY OF STATES 

In this section we briefly discuss the meaning and the implementation of the semiclassical approximation for a 
system of particles possessing one-body Hamiltonian H with the eigenvalue equation 

H\^n)=E^\^n). (4) 
where {t/jn) and En denote the n-th single-particle eigenstate and eigenvalue, respectively. 
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In view of studying the various thermodynamic properties of the system, one is usuaUy interested in evaluating 
sums of the form 

n 

where F is some function of the single-particle energy levels. The semiclassical approximation consists in considering 
the energy levels so close to one another that the sum over them can be taken as an integral over a phase space 
density: 

Y,nEn)= [9{E)F{E)dE. (6) 

To this end we define the density of states 

g{E)='£6{E-En), (7) 

n 

where S is the Dirac delta distribution. In the following we investigate for the parameters of the Table I whether it is 
justified to approximate the density of states semiclassically according to 

where i^(x, p) = p^/ (2M) + V(x) denotes the single-particle classical Hamiltonian. In systems where the exact energy 
spectrum is known, e.g. a harmonic trap, the validity of this approximation can be determined by comparing (7) and 
(8). However, in the case that the energy spectrum is unknown, as for the present anharmonic potential (2), we must 
take a diff"erent approach to justify the semiclassical approximation. To this end, we determine not only the leading 
semiclassical approximation go{E)^ but also its subleading correction gi{E). After setting the relevant energy scale 
of the system, the validity of the semiclassical approximation can be determined by the relative magnitude of the 
correction. We first treat the general case and, then, specialize to the anharmonic trap of the Paris experiment. 



A. Semiclassical Expansion 

To obtain the semiclassical expansion of the density of states of the anharmonic potential, we follow the general 
procedure worked out in Ref. [45], which derives the density of states in form of a gradient expansion 

g{E) = go{E) + g^{E) + ■ ■ ■ . (9) 

The first term corresponds to the semiclassical approximation ^o(^) iii (8), which yields after integration over the 
momenta 

The second term is the subleading correction to the density of states 

\ M " 

gi{E) = - 



2AM 



r(L>/2-2) 
xe[E- y(x)] VV(x) [E - F(x)]*"^ . 

(11) 



We note that, for the physically relevant dimension D = 3, the spatial integral (11) diverges due to the behavior of 

— —3 

the term [E — V{-x.)]'^ at the classical turning points. This problem can be solved through the method of analytic 
continuation as demonstrated later by applying (11) to the anharmonic potential (2). 

It remains to determine the energy scales for which the semiclassical approach is supposed to be valid. At high 
enough temperatures, the particles are classical. When the temperature is gradually lowered, or other parameters 
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are changed, classical statistical mechanics is no longer exact. At first, a semiclassical approach may work, but for 
certain conditions it may completely fail. The temperature, particle number, and trap potential parameter scales 
in which the semiclassical approximation is valid depend on the underlying statistics. For Bosons in the gas phase, 
the thermal distribution is characterized by the thermal energy /c^T. If this energy is sufficiently larger than the 
the energy of the single-particle ground-state, the leading semiclassical approximation is valid. At the BEC critical 
temperature Tc, however, the ground-state occupancy becomes macroscopic. This effect must be treated outside of 
semiclassics because it corresponds to a broken symmetry state associated to the divergence of the Bose distribution 
(e^^ — 1)"-^ for E = 0. We will see explicitly in the case of the anharmonic oscillator that the subleading correction 
is not negligible for energies near the ground-state energy. The thermal energy at the critical temperature is thus the 
relevant energy scale for determining the validity of this approach. 

For Fermions the considerations are slightly different. Due to the Pauli exclusion principle, as T ^ the particles 
in the system do not all drop to the ground-state but rather fill every energy level up to the Fermi energy Ef, which 
is defined as the energy of the highest occupied state. As the density of states g{E) typically increases with a power 
of E^ the majority of particles in the system will have energies much larger than the ground-state value and their 
average energy is of order Ep, even at T = 0, provided is large. Therefore, in principle, a spin-polarized Fermi gas 
can be well described at all temperatures in the semiclassical approximation provided the particle number is large 
enough. The relevant energy determining the validity of this assertion is simply the Fermi energy Ep. Numerical 
calculations indicate that this is indeed the case [46]. 



B. Leading Semiclassical Approximation 

To consider the anharmonic potential (2), we must first determine explicitly the leading and subleading terms of 
the semiclassical density of states. We can, then, apply the leading term to calculate the Fermi energy, and with this 
result return to verify the a priori assumption of the validity of the approximation. 

The form of the semiclassical density of states in the anharmonic potential (2) can be found from (10). Performing 
the spatial integral we obtain the result 



2(E + A^)i X^rjE X'^r]AE 

4{E^AE)i _ 
3e| 



^ > 



-AE <E<0 

(12) 



where it is meant in this and other equations, that the function takes on nonzero values for £' < only for super- 
critical rotations. The density of states (12) for the parameters of the anharmonic trap described in Table I is shown 
at several rotation frequencies Q in Fig. 2. For subcritical rotations, the density of states vanishes for < 0. For 
supercritical rotations, the density of states is nonzero also for —AE < E < 0^ due to the dip in the potential. 

There are two important limiting cases which shed some light on the behavior of the system. The first one is the 
harmonic limit ^ 0, for which we have from (12) 

which corresponds to the known result [5]. Note that the harmonic limit (13) is valid only for subcritical rotations; 
for supercritical rotations the density of states (13) becomes negative corresponding to an instability due to the fact 
that particles escape the trap. 

The second physically important case is the pure quart ic limit in radial confinement. It is realized by either r/ = 0, 
i.e. a critical rotation, or A^Ty, so that the density of states (12) takes the form 

g^^{E) = ^^, ri = 0,E>0. (14) 
3e| 

Examining (3) and (12), we see that the rotation parameter rj appears in the density of states only in ratios to the 
anharmonicity parameter k. For supercritical rotations, we will find that there is a symmetry between the range of 
rotations, rj — oo (Q oo) to 77 = (critical rotation), and the range of anharmonicities from small to ^ cxo. 
The supercritical properties of the system can thus be seen as an interpolation between the harmonic limit (13) and 
the pure quartic limit (14). 
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FIG. 2: (Color online) Density of states (12) with the parameters of Table I for the anisotropy and the anharmonicity of the 
trap. The solid lines are from bottom to top at rotations of Q = 0.1 o;^ to = 1.3 o;^ in steps of AQ = 0.2 lj±. The dashed 
line is the density of states at the critical rotation Q — uj± . 

C. Subleading Semiclassical Approximation 

We now calculate the subleading semiclassical approximation gi{E) as given in (11) for the particular case of the 
anharmonic potential. A direct integration in (11) for = 3 yields a divergent result. However, we can rewrite the 
spatial integral in terms of the Beta function 



P{x, y) = f 
Jo 



ds s""-^ {1 - sf-\ 



for which in Ref. [47] provides the identity 



/3{x, y) = 



r(x)r(t/) _ 

r{x + y) ' 



(15) 



(16) 



Thus, the formally divergent integrals can be treated as finite based on the analytic continuation of the Gamma 
function through the relationship T{u + 1) = uT{u). 

To evaluate the integral we need the Laplacian of the anharmonic potential, which reads: 



(17) 



With this and the full potential given in (2), we obtain from (11) the first correction to the semiclassical density of 
states: 



-1/2 



[2X^7] -1) {E^AE} 

{2X^7] -1) (E^AE)-^^^ 
24f^hl 



1 



E>0 



-A^ < ^ < 



(18) 



where, again, the lower line is valid only for supercritical rotations. 

To determine the effects of the trap parameters on the relevance of the correction term, we examine the ratio 
gi{E)/ go{E). We first see that it is in general reduced in both the sub- and supercritical regimes as {E + AE)/ez 
increases. Likewise, the importance of the correction term decreases as the anharmonicity decreases. As the 
correction is largest for low energies, the most interesting case for supercritical rotations is when E < 0^ where we 
have 



gi{E) 2X^r]-l 



go{E) 



32 



E^AE 



-AE <E <0. 



(19) 
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FIG. 3: (Color online) The semiclassical density of states (12) and the leading semiclassical correction (18) for several rotation 
frequencies, using the parameters of Table L 

From (19) we see that for a fixed E-\-AE^ the relevance of the correction increases linearly with the rotation parameter 
rj. Furthermore, we see that the correction (18) diverges as £^ ^ —AE. However, this divergence is only of the form 
{E + A£^)~^/^, and so any interesting integral of the form (6) over the density of states will still have a finite 
contribution from this correction, provided F{E) is not itself divergent. For Bose Systems, F{E) contains a factor 
of [exp/3(/i — E) — 1]~^. The phase transition to Bose-Einstein condensation happens as /i ^ £^0, where Eq is the 
lowest energy eigenvalue. Thus, in the condensate phase, an integral of the form (6) has an improper lower limit for 
F(yE) blows up as £^ ^ —AE^ since it is the lowest available energy. This problem can be overcome by dividing the 
gas explicitly into a condensate plus thermal particles. 

We can now check the validity of the semiclassical approximation in the case of the anharmonic trap. Both the 
leading (12) and the subleading (18) term are plotted in Fig. 3 for the parameters of Table I , where we see that even 
for supercritical rotations the subleading correction becomes negligible on energy scales of more than a few Cz. In 
the following section we determine the Fermi energy Ep of the system, from which we will see that, for the chosen 
parameters, the difference between Ef and the minimum potential is orders of magnitude larger than e^. With this, 
we justify the semiclassical approximation (12) to the density of states. 

IV. GROUND-STATE PROPERTIES 

Though not fully accessible from the experimental point of view, the ground-state and its properties are crucial 
for understanding the system. The Fermi energy, e.g., carries information about aspects like dimensionality, trap 
geometry, and particle number and, thus, providing information necessary to handle important physical issues. 

In this section, we consider some properties of the ground-state like the dependence of the Fermi energy on the 
particle number, the ground-state energy, the spatial and momentum distribution and, lastly, the ballistic expansion 
of the system. 

A. Fermi Energy 

The Fermi distribution is given by 

f{E, /3,M)= ^,(J,)^^ , (20) 

where /i is the chemical potential and /3 = l/ksT with ks the Boltzmann constant and T the temperature of the 
system. At T = 0, i.e., ^ oo, the Fermi distribution becomes a step function at the Fermi energy, which is defined 
as the chemical potential at T = 0: 

lim f{E, (3, Ai) = eiEp - E); Ep = Ai|r=o • 
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The ground-state particle number is therefore given by 



N{Ef) = j " 9{E)dE. (21) 
Evaluating this integral with the semiclassical density of states from (12) yields 



N{Ef) = 



4(^i. + A^)i \^r]\El AEEf 2AE 



15e| 



4e2 3e2 ISe? 



(22) 



for the first case, in which one has Ep > 0. For the second case, in which —AE < Ep < 0, the particle number is 
given by 

15e| 

The inverse of (22-23) determines the Fermi energy as a function of the trap parameters and the particle number. In 
the harmonic [5] and pure quartic limits we have 

2 

= {QNX^^^ e,, eQ"= (^K^iv)%„ (24) 

respectively. We note that the Fermi temperature E^^ /ks is of the same order of magnitude as the critical temper- 
ature for non-interacting Bose-Einstein condensates in the quartic regime [14] 



QU 



2K.2N 
C(5/2)^^ 



(25) 



with C(5/2) ^ 1.342. 

It should be pointed out that negative Fermi energies are allowed by (22-23) for supercritical rotations, in which 
case the gas is confined entirely to the dip in the potential, creating a hole in the center of the trap analogous to the 
one found in two-dimensional BECs [48]. In this interesting regime (— < Ep < 0) we see that the Fermi energy 
also takes on a simple analytic form, 

2 

E^"" = -AE + 6, (^^'^n] ' , (26) 



where the superscript DO stands for donut, in reference to the shape that the gas assumes when all the particles are 
in this energy interval, as discussed below. Notice that in the 'donut' regime, for a fixed particle number the Fermi 
energy is set a distance above the trap minimum that is independent of the rotation frequency, but is directly related 
to the anharmonicity and to the number of particles. 

We now have three cases where the Fermi energy takes a simple analytic form corresponding to the three interesting 
rotation regimes: the harmonic limit for subcritical rotations, the pure quartic limit for critical rotations, and the 
'donut' regime for supercritical rotations. Although we can not obtain from (22-23) a general explicit expression for 
the Fermi energy, we can determine its behavior from these special cases. In Fig. 4 we see that for low Ep, the general 
supercritical behavior approaches the harmonic oscillator behavior. As rotation increases to the critical value, the 
curves approach the pure quartic limit. As the rotation becomes supercritical, the curves are offset to more negative 
energies and the 'donut' regime eventually sets up. 

In Fig. 4 we also see that for any Fermi energy possible, the quantum correction term to the semiclassical approxi- 
mation (18) is negligible, thus justifying our a priori use of the semiclassical approximation. 

In the following we use the Fermi energy at zero rotation as a convenient energy scale. For our parameters it is 

Efo ^ 430e^. 

The Fermi temperature Tp = Ep/ks at zero rotation provides a corresponding temperature scale, 

Tfo = Ep/kB - 230 nK. (27) 
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FIG. 4: (Color online) Relation between Fermi energy and particle number given by (22-23). The solid curves are for subcritical 
rotations Q/cj± = 0, 1/3, 2/3 from top to bottom on the left side, and for supercritical rotations Q^/uj^ = 1.1, 1.3, 1.5 from top 
to bottom on the right side. The dashed curve on the left is the harmonic limit (24) at zero rotation. The dashed-dotted curve 
on both sides is the pure quartic/critical rotation limit (24). The Fermi energies in the 'donut' regime are given by (26) and 
are bold-faced. The parameters of the Table I were used, and a guide-line dX N — ?>x 10^ indicates the chosen particle number. 



B. Internal Energy 



The total internal energy of the system in the ground-state is given by 

U{Ef) = [ ' Eg{E)dE. 

J — OO 

Evaluating this integral with the density of states from (12) yields 

W{Ef) X^V 



U{Ef) 



E% ^ AEEj. 8AE^ 

6^ 356^ 



e2 

L 



(28) 



(29) 



2 

W{Ef) 

with the remark that the upper hne is vahd for Ep > and the lower one for —AE < Ep < and the abbreviation 



W{Ef 



8Ef (Ef + AE) 5 16 (Ef + AE) ^ 



15e| /t2 



105e|K2 



(30) 



In the harmonic and purely quart ic limits, respectively, the internal energy it is related to the particle number and 
Fermi energy by 



^HO ^ 3^^HO ^QU ^ ^^^Qu 

4 ^ 7 ^ 



There is also a simple relation in the 'donut' regime, 



U 



DO 



7^- 



DO 



7 



-AE N. 



(31) 



(32) 



We see in Fig. 5 that the internal energy (29) for subcritical rotations is bounded by the ^ harmonic oscillator 
limit and the ^ oo pure quartic limit (31). This is due to the symmetry between the anharmonicity and the 
rotation, which causes the latter to enhance the former. 



C. Particle Density 



To obtain the particle density and the momentum distribution in the semiclassical approximation, we start from 
the Wigner function [49] in the semiclassical approximation 

^(^'P)=l + exp{/J[iJ(x,p)-,]}' 
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FIG. 5: (Color online) Ratio of internal energy (29) per particle to Ef + for fixed N and the parameters of Table I. The 
dashed horizontal line is the harmonic limit and the dot-dashed horizontal line is the quartic limit (31). The vertical lines 
indicate the critical rotation (Q = (jJ±) and the start of the 'donut' regime (32), where Ep < 0. 



which delivers the particle density through 



and the momentum distribution through 



.(x) = / 



n(p) 



z/(x, p) 



^(x, p). 



(27rn)3 

At T = the Wigner function becomes a Heaviside theta function and is given by 

i.(x,p) = e[Ef -F(x,p)]. 



(34) 



(35) 



(36) 



In this limit, both the particle density and the momentum distribution assume a well-defined shape. The particle 
density is given for a general velocity independent potential T^(x) by 



n(x) 



3 



e [Ef - T^(x)] . 



(37) 



Prom (37), we see that a given rotation $7do> at which the hole in the trap center appears, the system looks like a 
'donut', so to speak. The value of f^DO can, therefore, be determined from (26) by setting Ep = and yields 



DO 



A2 



(38) 



In an effectively two-dimensional (uniform over an axial length Z), interacting BEC trapped by the radial components 
of the anharmonic potential (2), rotation can lead as well to a hole in the center of the trap [48], where, in the Thomas- 
Fermi approximation, a similar analytic expression can be written down: 



(l2n^Na/Z)l 



A2 



(39) 



where a denotes the s-wave scattering length. In the case of non-interacting bosons, Vt\^ coincides with uo^^ but this 
is not true for fermions, where an effective interaction arises from the Fermi statistics and I^do does not equal uj±_ 
even in the absence of interactions. A discussion of this Pauli "pseudopotential" in the case of harmonically trapped 
Fermi gases is provided by Ref. [41]. 

The value of I^do could also be obtained from the length scales Bi^^ set up by the Fermi energy, describing the 
extent of the gas in space. 



[y(x)-i?F]U=flTP = o, 



(40) 
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where i, j range over z. These length scales describe the maximum extent of the coordinate space distributions 
at T = and are called the Thomas-Fermi radii Rj^ . For the anharmonic potential (2), they are given by 



= ^^jEF+ei-rj)AE], (41) 

where the positive root Rj^ gives the outer radius of the distribution and the negative root has a physical meaning 
only in the 'donut' regime, where it gives the inner radius. We note that for supercritical rotations {Ef < 0), the 
density distribution has the shape of a shell, of which we then designate the outermost limit in the z = 0-plane 
as the Thomas-Fermi radius Rj^ and the largest extention in the 2:-direction as Rj^ . In the harmonic limit, the 
Thomas-Fermi radii are 

In particular, we are interested in the optical density which can be directly measured by absorption images. Per- 
forming an integration over the axial ^^-direction, we obtain 



n{r±) = J 



2 



"K.^M-^^^^S^. (43) 



The rotation dependence of the T = column density (43) is shown in Fig. 6 for the anharmonic potential (2). Just 
above the critical rotation the distribution becomes nearly homogeneous, and at high supercritical rotations the gas 
eventually enters the 'donut' regime, where the gas is entirely confined to the dip in the trap. In contrast, without the 
quart ic trapping the gas spreads out to become an infinite pancake as the rotation approaches the critical rotation, 
as we can see from the divergence of the harmonic oscillator Thomas-Fermi radius (42) at the critical rotation. 



D. Momentum Distribution 



The Fermi energy also sets a momentum scale through the Fermi momentum pr = \/'lM \Ef + Q[—r])^E\. It 
represents the maximum extent of the momentum distribution at T = 0. For the anharmonic trap we have azPF = 

A simple dependence on the particle number emerges again in the harmonic limit 



Vf 



(44) 



At T = the momentum space distribution (35) for the anharmonic potential (2) reads 



n(p) 



^4 
167r2fi3 



(g + l)tan-i(^)- V>0 
(5 + 1) [7r-tan-i(^)] q>0 
{q + l)7r; q<0 



with the dimensionless variable 



AE 



2MEf 



(45) 



We see that, in contrast to BEC, where the momentum distribution of the condensed atoms inherits the anisotro- 
py of the trap, the fermionic momentum distribution of a non-interacting gas remains isotropic despite the trap 
anharmonicity. This feature is drastically affected by presence of the anisotropic dipole-dipole interaction, as then 
the Fock term of the two-body mean-field energy induces a characteristic deformation of the momentum distribution 
[50]. Note that we have only considered the co-rotating frame; as we will see, the rotation leads to an anisotropy in 
the momentum distribution in the lab frame. 
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FIG. 6: (Color Online) The T = spatial density distribution (43) for the parameters of the Paris experiment at increasing 
rotations. The anharmonic term prevents the formation of a 'pancake' at near-critical rotations and causes the gas to take on 
a 'donut' shape for large supercritical rotations. 



E. Free Expansion 

Measurements in cold gases are usually carried out by turning off the trap and shining a resonant laser on the 
sample after some time of free expansion. These are called time-of- flight measurements (TOF). 

So far we have treated the problem in the corotating frame. In considering TOF measurements, it is necessary 
to consider the laboratory frame, where the gas expands ballistically. The transformation of the momenta from the 
corotating frame to the lab frame yields 

Piab = P + Mnr4>, (46) 

where ^ is the tangential unit vector. Since the trap is radially symmetric, the coordinate transformation of the 
position does not affect the distribution. We have then 

^lab(xiab, Plab) = ^(x, p) = niab(xiab, Plab - MVtr^). (47) 

Thus although the momentum distribution does not inherit the anisotropy of the trap, the trap rotation leads nonethe- 
less to an anisotropic momentum distribution. 

To study the expansion in detail we apply the ballistic substitution 

(x,p)^(x-H^,p), (48) 
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which yields a relation between the lab frame distribution and the corotating distribution: 



^lab(xiab,Plab,^) = ^ 



Xlab 



M 



,Piab-Mnri"Vo 



(49) 



where r\_ and 4>o denote respectively the radius and tangential vector at the position xq = xiab — , and can be 
expressed as 



^^^00 = Z X Xo = Z X I Xlab 



Plab^ 

M 



(50) 



with z the unit vector in the axial direction. 

Of particular experimental interest is the density distribution at a time t after the trap is turned off, 



2M 



-V X 



Plab A 

M ) 



(51) 



For the case of harmonic trapping and zero rotation, the kinetic term can be combined with the harmonic term and 
the integral over momentum space in (51) takes nearly the same form as the integration over momentum space in 
(8). Unfortunately, in the case of the anharmonic potential (2) with rotation, the integration can not be cast into a 
simple form and does not provide an insightful analytic expression. 
A measure of the spatial density is given by the aspect ratio 



e{t) 



with the respective widths 



Wi{t) 



?^lab(x, t) 



(52) 



(53) 



For harmonically confined Fermi gases with zero rotation, it is known that anisotropic samples become asymptotically 
isotropic for large expansion time [51], since their widths scale in time as ^^^^{l + (jj^ft^)^ and the aspect ratio as 



e{t) 




1+6^2^2- 



(54) 



In the case of a rotating gas with additional quartic trap, it is interesting to see how the rotation affects the 
asymptotic value of the aspect ratio, which reflects the lab frame momentum distribution. In Fig. 7 the aspect ratio 
(52) has been evaluated numerically for the anharmonic potential. Comparing to the harmonic trap at zero rotation, 
we see that the anisotropy of the momentum distribution increases with the rotation frequency for the anharmonic 
trap. It is worth mentioning that the value of e{t) for large asymptotically large t depends on the rotation, so that 
the aspect ratio could be used to measure the actual rotation frequency in an experiment. 



V. TEMPERATURE EFFECTS 



After having studied the effects of rotation on the ground state of an anharmonically trapped polarized Fermi gas, 
we turn our attention to the thermodynamic properties of such a system. The standard thermo dynamical recipe is 
specialized to the case of the anharmonic potential and the consequences for different rotation frequencies are explored 
in both low- and high-temperature regimes. 



A. Standard Grand- Canonical Approach 



In this section we describe the thermodynamical properties of the system within the framework of the grand- 
canonical ensemble, where the most fundamental quantity is the grand- canonical partition function. For particles 



14 




FIG. 7: (Color online) The aspect ratio (52) of the spatial density for the parameters of the Paris experiment at different 
rotation frequencies as a function of expansion time. The dashed curve is for the harmonic trap at zero rotation, and the solid 
curves are for the anharmonic trap at rotations of Q/o;^ = 0, 0.75, 0.99, 1.25 from top to bottom. 



obeying Fermi-Dirac statistics it is given by 



1 + e- 



(55) 



where En, as before, is the eigenvalue corresponding to the eigenstate |n). 

The grand-canonical free energy is obtained from the grand-canonical partition function through 



(56) 



Inserting (55) in (56) and making use of the semiclassical approximation, i.e., writing the discrete sum as an integral 
over the density of states (7), we obtain 



g{E)\n 



1 + e-''^^-^) 



dE. 



(57) 



Performing a partial integration of (57), one sees that the integral for the grand-canonical free energy contains the 
Fermi distribution: 



G{E)- 



dE 



where G{E) denotes the number of single-particle states with energy less than E 

nE 



G{E) = / g{E')dE'. 

J — oo 



(58) 



(59) 



In view of using the semiclassical approximation for the density of states in (12), we define the incomplete Fermi 
integral 



1 + e^z' 



(60) 



which is called incomplete due to the fact that its integration range is smaller than the one of the original Fermi 
integral, defined by 



(61) 
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For convenience we introduce the further shorthand 



(62) 



With (12) and (58), we write the grand-canonical free energy of a system contained in the anharmonic potential 
(2) in terms of the following Fermi integrals: 



where a{r]) is the sign of 77, defined as 



+ ,inc 
7 
2 



\6 3/-+4nc A 10 5/-+,inc 



air]) 



320/3^3 



+1; ri>0 
-1; v<0 



(63) 



(64) 



To obtain the other thermodynamical quantities, we express the grand-canonical free energy in terms of a Legendre 
transformation 



with the differential 



dTc = dU- SdT - Ndfi. 



(65) 



(66) 



Therefore, the particle number N and entropy S can thus be determined from the grand- canonical free energy 
according to 



TV: 



dJ^G 



dT 



(67) 



For the particle number we obtain 



TV 



2e(-r7)a + a(77)C 



+ ,inc 



2/35e|/t2 
24/3e^K2 ' 



2/32e2K 



(68) 



and for the entropy we get 



TS 



5V^e(-r?)Ct 



+ - 



+ ,inc 



2/35 el K 2 



The internal energy can then be obtained from (65) yielding 

5V^ [2e(-77)C^ + ^('7)C^'"'' 



+ 



A^r?3c+-'- 
6/32e, 



4/32 e|K2 
NAE. 



6/32e, 



(69) 



(70) 
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FIG. 8: (Color Online) Heat capacity per particle at constant particle number and trap parameters for varying temperatures 
and rotation speeds from (73). The energy scale is Efo, the Fermi energy of the system at zero rotation = 0). The 
temperature scale is the Fermi temperature, T/o = Efo/kB- The angular frequency scale is the critical frequency (jJ±. 



Notice that, in the pure quartic limit (ry = 0), a simple relation between the free and the internal energies arises 



The heat capacity of the system at constant particle number and trap parameters is given by 

dU 



C 



dT 



(71) 



(72) 



N 



To evaluate this derivative, we must, as usual, obtain the temperature derivative of the chemical potential from the 
With (70) we then obtain 



35V^ 2e(-77)ct + cT{'n)^Y'^ 



2e(-r7)C+ + ^(r?)ct'^'' 



f32 n2 



3^2^^+, inc ^6^3^+, inc 



2e(-r?)CT + ^(r?)C' 



/3 

+ ,inc 



5^1 



/32e2, 



(73) 



Since this analytic expression is not particularly illuminating, we analyze the behavior of (73) for a range of 
temperatures from Fig. 8. At low temperatures the heat capacity is nearly linear in temperature, and at high 
temperatures the heat capacity approaches C = ^Nksf^ for a fixed rotation. At high temperatures a dip appears 
with a minimum occurring at a rotation that increases with the temperature. This interesting case and the general 
high-temperature behavior of the heat capacity will be discussed in detail below. 



B. Low- Temperature Approximation 



In the low-temperature regime, we can understand temperature effects by using the Sommerfeld expansion as 
sketched in Ref. [54] . To this end we consider the incomplete Fermi function with X = /i + 

ou poo 

C+'"^(e^^,/?Ai?) = -^ / E''-'fiE,p,X)dE, (74) 
i ['^) Jae 
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where f{E, /3, X) is the Fermi distribution (20). A partial integration yields 



Tiu) 



1 



f{AE, (3, X)AE- 



pOO 

/ E 

Jae 



where the derivative of the Fermi distribution is given by 

df{E, P, X) 



,df {E, 13, X) 
dE 



dE 



dE 



(75) 



(76) 



Consider the case that \E — X\ ksT. We see that the derivative (76) becomes very smaU. If X > 0, then the 
lower hmit of integration in (75) can be extended to —oo and f{AE^ X) ^ 1. If X < then the integrand vanishes 
and f{AE, X) ^ 0. We now expand E^ in a Taylor series about E = X, 



E"" 



and integrate (74) with (77) to obtain 



E 

n=0 



r(^ + i) 

V{v+l-n) 



X' 



{E - XY 
r(n + l) 



-2n 

"2n) 



(77) 



(78) 



It is worth remarking that the Sommerfeld expansion (78) is valide for X > and vanishes for X < 0. Furthermore, 
the an are constants given in terms of the Riemann zeta function C,(n) as 



1 



C(2n). 



22(n-l) 

For the case of the complete Fermi integral, Eq. (78) reduces to the usual Sommerfeld expansion given by 



C(e''^) 



oo 



an {(3X) 



-2n 



^ r(z/ + 1 - 2n) ' 

n=l ^ ^ 



(79) 



(80) 



for X >0. 

Although the series in (78) and (80) terminate for v = 0,1,2..., the expansion is not exact since the limits of 
integration were extended in (74) . 

It is useful to obtain an explicit equation for the temperature dependence of the chemical potential using the 
Sommerfeld expansion (78) and the particle number equation (68). To this end, we substitute /i = Ep + 6 and keep 
terms linear in 6 as well as quadratic in temperature to obtain up to third order contributions 



Ep 



SQgoiEp) 



6(Sf + AS)5 3X^n 



(81) 



where ^o(^f) is the density of states (12) evaluated at the Fermi energy. We see that for low temperatures the 
chemical potential decreases with the square of the temperature. This result is compared to the numerical inversion 
of the particle number equation (68) with both the exact evaluation of the Fermi integrals and with the evaluation 
using terms quartic in temperature from the Sommerfeld expansion in Fig. 9. 



C. High- Temperature Approximation 

In the high-temperature limit we expect quantum statistics for both bosons and fermions to approach classical 
statistics. Here we derive the high-temperature expansions of the thermodynamic properties of a Fermi gas in the 
anharmonic trap of the Paris experiment. We also consider in more detail the unique rotation dependence of the heat 
capacity in the high-temperature limit. 
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FIG. 9: (Color online) Evaluation of chemical potential from particle number equation. A numeric inversion of (68) is shown 
for the general case. For the low temperature expansion (81) is shown, and for the high temperature expansion (87) has been 
solved analytically for /a. 



A high-temperature expansion of the incomplete Fermi integrals is obtained as follows. We rewrite (60) as 



c 



e -^z 
1 + e~^z 



-dx, 



(82) 



where z = e^^^^^^\ As seen in Fig. 9, if the temperature is high enough, e^^ < 1 and, thus, e (^^^ z < 1 so that we 
can use the binomial expansion for the factor (1 + e~^z) ^ . Then we obtain 



C+''"=(^,/3A£:) = ^Q(j.,i/3A£:) 



with the definition 



Q{v,l3AE) 



^(^ 



1 

(^) Jf3AE 



(83) 



(84) 



where T{x) is the Gamma function. 

For z ^ in the high-temperature limit, it is sufficient to keep only the ffist term in the sum, yielding 

lim C^''"' = lim C^''"'(e^^^+^^\/5A^) 
For positive integer n, Q{n,j3A.E) can be expanded exactly: 



(85) 



(86) 



The high-temperature expansion is compared to the exact evaluation of the Fermi integrals in Fig. 9 in the context 
of a numerical inversion of the particle number equation to solve for the chemical potential. The approximation is in 
general very accurate for temperatures larger than half the Fermi temperature. 
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From (85) we can determine the high-temperature behavior of the general expressions for the thermodynamic 
properties of the system. For the particle number, we get 

^Pi^+AE) r 

2e(-J7) + <t{v)Q 



N 



+■ 



(87) 



It should be remarked that (87) is to be understood as giving the temperature dependence of the chemical potential 
in the high-temperature regime, so that /i is changed to maintain a constant value of N on the left-hand side of the 
equation. Note that the high-temperature chemical potential can be determined explicitly by taking the logarithm of 
(87). 

Now we are in a position to write down the high-temperature versions of the internal energy (70) 



U 



NknT 



and of the heat capacity (73) 



C 



■(3AE + 



[2e{-v) + a{v)QilfiAE)] 



+ 



(88) 



-1/2 -/3A_E 



[2e{-r,)+a{r,)Q{lpAE)] 



n[2e{-r]) + a{rj)Q{lpAE)y 



+ ■ 



(89) 



The high-temperature behavior of the heat capacity (89) is shown in Fig. 10. For subcritical rotations we have 
C > 2.5 NkB, for critical rotations C = 2.5A^/cb, and for supercritical rotations C < 2.5 Nks- For a fixed rotation, 
the heat capacity approaches in all cases C = 2.5 Nks as T ^ oo. For a large constant temperature, the heat capacity 
also approaches C = 2.5 Nks in the high rotation limit, Q/uj_\_ oc. 

The most interesting feature of the high-temperature heat capacity is the dip found in the supercritical region. 
Differentiation of (89) with respect to Q shows that the minimum heat capacity occurs at fSAEco ~ 0.948. Inter- 
estingly, the value of the minimum heat capacity, Cq ^ 2.31SNkB^ is independent of temperature and rotation. The 
relation between the temperature of the system and the rotation Oq corresponding to the minimum heat capacity can 
be written more clearly as 



■ ^_L1 



(90) 



The rotation Qq at which the minimum heat capacity is located increases monotonically with temperature and has 
the asymptotic form fto oc T^^^. Thus, although the minimum is present at any temperature, it occurs at increasingly 
high rotations so that the heat capacity flattens to C = 2.5 Nks as T ^ oc for any finite rotation. Furthermore, the 
similarity of (90) to (38) and (39), which give the critical rotation for the transition to the zero temperature donut 
regime, suggests that the minimum of the heat capacity is related to the transition of the gas into a configuration 
where the particles are mostly contained in the minimum of the supercritical potential. Recalling that AE is the 
potential height at the trap mininum, one should expect that the "donut regime" is realized when the thermal energy 
of the particles is of the order of AE. Thus, setting ^ ksT should provide an estimate of VLq. Indeed, this gives 



(91) 



which is not bad at all. 



D. Virial-Theorem Approach 

At high enough temperatures, the quantum mechanical effects may be neglected. Thus, the high-temperature 
properties of systems in an anharmonic trap can also be explained in the context of classical statistical mechanics. To 
do so, one has to average the desired quantities over an ensemble. For a function 0(x, p) of the phase space points 
such an average {0(x, p)) is given explicitly through 

_ /d3.d3pQ(x,p)e-/^^(-^P) 
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FIG. 10: The high-temperature behavior of the heat capacity as a function of rotation, evaluated using the parameters of the 
Table I and T = 1 /xK with (89). From Eq. (90), one obtains that the minimum of the heat capacity at T = 1 /iK is located at 
Qo ~ 3.43 ix}±. 
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Cv = 2.318NkB' 






pAEco = 0.948 



where (x, p) is again the classical Hamiltonian. 

In this context, the classical virial theorem states that for N particles in thermal equilibrium at temperature T in 
a potential the classical ensemble averaged potential energy (V) and translational kinetic energy (K) of a particle 
in the system are related according to [53]: 

1 / dV \ 1 

(i^.) = ^(^-9=2^-^' (93) 

where Xi are the components of the position vector of a particle and K = Ki. For a cylinder symmetrical potential 
of the form ]/ = V(r, z), one has 



The anharmonic potential, as given in (2), can be decomposed according to 

V{r^,z) = VrAr±) + V.{z) (95) 

with the terms 



VrAr±) = I 



2„' ± 



VM = ^ (96) 



2 

Using (94), we obtain for the z-degrees of freedom the familiar harmonic potential result 

(i^z) = {V.) ■ (97) 
For the radial degrees of freedom we find correspondingly 



(/f.) + {Ky) = 2 {V;j - A^^l^. (98) 



The total internal energy of a system of N non-interacting particles in the potential V is then given by 

V = N{{K) + {V)) = ^NkBT + ^l^^'fjl^ (99) 
2 4 

Before we consider the general result, let us make some remarks about Eq. (99). By using the radial symmetry of 
the problem, one can express the average square radius through 
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with the momenta 



rui 



f 

Jo 



p^dpe 



(101) 



With the help of Ref. [47], the hmiting cases can be analyticaUy explored. In the case of a pure quartic trap (77 = 0), 
the average radius is given by 



showing, as expected, that tv sets up the important length scale. 
For the case of a pure harmonic trap, i. e. = 0, one has 



2 \HO 



{rl) 



(3Muj\r] 



(102) 



(103) 



From (103), we conclude that the average position of the particles diverges as 77 approaches 0. Nevertheless, the 
energy per particle, as given from (99), remains SksT. 
In general we obtain 



(rl) 



2at 



(K/3e,) 



r{-cT{rj){/3AE)i- 



-PAE 



[20{-ri)+a{7i)Q{l,pAE)]^ 



(104) 



Substituting this into (99) we recover (88). We, thus, see that the classical statistical treatment of the system is 
equivalent to the high-temperature expansion which we have developed by keeping only the leading term in (83). 



VI. CONCLUSIONS AND OUTLOOK 



We have considered non-interacting spin-polarized Fermi particles in an anharmonic trap in the presence of a 
centrifugal potential due to rotation of the system. Applying a semiclassical treatment to the density of states, we 
have worked out the temperature dependence of several t her mo dynamical properties. At absolute zero we have studied 
the density of particles and the momentum distribution for an anharmonically trapped system. The analysis of the 
particle density shows that a well defined geometry characterizes the ground state and has revealed a tendency for the 
particles to move away from the center. An analytical expression was obtained for the value of the rotation frequency 
Odo at which a hole is created at the center of the trap, giving rise to a regime where the gas is 'donut '-shaped. The 
frequency ^do was shown to depart from uj±, which characterizes the stability limit for harmonic confinements, due 
to the interplay between the fermionic statistical interaction and the anharmonicity, in close analogy to repulsively 
interacting bosons. At finite temperature, expressions for the chemical potential, the internal energy, and the heat 
capacity were given. These were followed by discussions of their dependence on both anharmonicity and rotation 
frequency. At high- temperature, a virial theorem approach was presented and compared to previous results as well 
as shown to be particularly useful in the intuitive limits of pure harmonic and quartic radial confinements. 

The most natural extension of this work would be the inclusion of interactions. In the low-temperature regime, 
short-range interactions are strongly suppressed for fully polarized fermions. To overcome this problem, one could 
include other Zeeman states by using an optical trap [55], yielding a spinor Fermi gas. Issues as (anti-)ferromagnetism 
and magnetization dependence of the different quantities [37] could be addressed in the presence of the additional 
quartic confinement. As for long-range interactions, no further Zeeman states are necessarily required. Due to its 
anisotropy, the dipole-dipole interaction is particularly interesting. Atoms with non-zero nuclear spin are possible 
constituents for systems with large magnetic dipole moment m. Degenerate stable, for ^^Cr (m = Qps) [56], and 
meta-stable states, for ^''^Yb (m = Sps) [57], are believed to be achievable. To date, only studies in two-dimensional 
systems consisting of a few atoms have already been performed [58], so that a three-dimensional many-body analog 
could be thought of. 

A further, and even simpler extension of the present work could be achieved by investigating the two-dimensional 
limit {ujz of the present system. By an appropriate tuning of the quartic constant, a quantum Hall regime 

could be realized making possible the observation of underlying Landau-levels through a plateau structure in the 
particle density when Q ^ u;_\_ [38]. For higher rotation frequencies, one would expect this structure to disappear and 
a hole to be created at a frequency value analogous to ^Ido- 
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